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The saturation length of aeolian sand transport ( L s ), 
characterizing the distance needed by wind-blown sand to 
adapt to changes in the wind shear, is essential for accu¬ 
rate modeling of the morphodynamics of Earth’s sandy land¬ 
scapes and for explaining the formation and shape of sand 
dunes. In the last decade, it has become a widely-accepted 
hypothesis that L a is proportional to the characteristic dis¬ 
tance needed by transported particles to reach the wind 
speed (the “drag length”). Here we challenge this hypoth¬ 
esis. From extensive numerical Discrete Element Method 
simulations, we find that, for medium and strong winds, 
L a oc V a /g, where V s is the saturated value of the average 
speed of sand particles traveling above the surface and g the 
gravitational constant. We show that this proportionality is 
consistent with a recent analytical model, in which the drag 
length is just one of four similarly important length scales 
relevant for sand transport saturation. 

1. Introduction 

Aeolian transport of sand occurs when a sufficiently 
strong wind blows over a sand bed [Bagnold, 1941; Shao, 
2008; Duran et al., 2011; Kok et al., 2012]. The two dom¬ 
inant transport modes are saltation, referring to particles 
hopping along the sand surface in characteristic trajecto¬ 
ries [Bagnold, 1941], and creep, referring to particles rolling 
and sliding along the sand surface [Bagnold, 1937]. Wind¬ 
blown, initially flat sand beds may evolve into bedforms, 
such as ripples and dunes, due to different kinds of insta¬ 
bilities [Claudin and Andreotti, 2006; Andreotti et al., 2010; 
Fourriere et al., 2010; Parteli et al, 2011; Charru et al, 
2013; Duran et al, 2014a]. 

For instance, dunes are thought to form due to an aerody¬ 
namic instability, namely a slight phase difference between 
topography and wind shear maxima on a periodically per¬ 
turbed sand bed [Jackson and Hunt, 1975; Hunt et al., 1988; 
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Kroy et al., 2002a, b]. If this phase difference is larger than 
the phase difference between sand transport and wind shear 
maxima, which is characterized by the saturation length 
(L a ) [Sauermann et al., 2001; Parteli and Herrmann, 2007a; 
Andreotti et al, 2010; Pahtz et al., 2013, 2014], the pertur¬ 
bations grow. This is one of the reasons why L s plays an im¬ 
portant role in aeolian dune formation. Indeed, L s controls 
the length of the smallest (“elementary”) dunes evolving 
from a flat sand bed and the minimal size of crescent-shaped 
barchans [Parteli et al., 2007; Claudin and Andreotti, 2006; 
Fourriere et al., 2010]. In contrast, the steady state dune di¬ 
mensions are controlled by the aerodynamic roughness ( z *) 
[Pelletier, 2009]. L a is also a key parameter in morphody- 
namic models of Earth’s sandy landscapes, such as aeolian 
dune models [Kroy et al., 2002a, b; Schwdmmle and Her¬ 
rmann, 2003; Parteli and Herrmann, 2007b; Narteau et al., 
2009; Parteli et al, 2009, 2014]. 

It has been a challenging task to predict L s as a function 
of wind and particle parameters, such as the wind shear 
velocity (u„), the mean particle diameter (d), the particle 
(p p ) and fluid density (p/), and the kinematic air viscosity 
(u). In fact, the main difficulty has been to understand 
which of the involved relaxation mechanism is the slow¬ 
est and thus the most important one. Sauermann et al. 
[2001] derived an expression for L s based on the assumption 
that L s corresponds to the length needed to eject particles 
from the sand bed. Based on measurements of the size of 
both subaqueous and aeolian barchan dunes, Hersen et al. 
[2002] proposed that L s is proportional to the drag length 
(Ld = ( Pp/pf)d), which characterizes the distance trans¬ 
ported particles need to reach the flow speed. Estimations 
of the wavelength of elementary dunes by Claudin and An¬ 
dreotti [2006] and measurements by Andreotti et al. [2010] 
later supported this proposition. Indeed, these measure¬ 
ments confirmed that L s is approximately proportional to d 
and essentially independent of it», as predicted by L s oc Ld- 
However, there is considerable room for alternative inter¬ 
pretations of these measurements. The analytical model for 
the saturation length of both subaqueous and aeolian parti¬ 
cle transport recently proposed by Pahtz et al. [2013, 2014] 
is also consistent with these measurements (without fitting), 
even though the model predicts that L s varies with u*. In 
this model, the four potentially most important relaxation 
mechanisms are all accounted for (ejection of bed particles 
and particle deceleration in particle-bed collisions, fluid drag 
acceleration of particles, relaxation of the fluid speed), and 
it turns out that neglecting any of them entirely changes 
the model predictions [Pahtz et al., 2014]. This shows that 
the identity of the most important relaxation mechanisms 
remains an open problem. 

Here we use Discrete Element Method (DEM) simulations 
for the particle phase to investigate aeolian sand transport 
saturation. This modeling technique considers interparticle 
interactions above and also with the sand bed and is thus 
more realistic than older modeling techniques [e.g., Almeida 
et al., 2007, 2008; Kok and Renno, 2009], which usually con¬ 
sider the sand bed as a flat, rough wall. However, it is also 
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computationally more costly, which is the main reason why 
this technique had not been used for modeling particle-laden 
flows until a few years ago [Carneiro et al, 2011; Duran 
et al, 2012; Carneiro et al., 2013; Duran et al., 2014b, a; 
Schmeeckle, 2014; Pahtz et al, 2015]. From our simulations, 
we find that the total mass of particles transported above 
the sand bed ( M ) relaxes significantly slower towards its sat¬ 
urated value (M s ) than the average particle velocity above 
the sand bed (V = Q/M, where Q is the sand transport 
rate above the sand bed) towards its saturated value (14), 
indicating that the drag length is not the dominant satura¬ 
tion length scale. Moreover, we find that L s oc Uj 2 /g, where 
g is the gravitational constant, when u» > 4-ut, where Ut is 
the dynamic threshold of sand transport (i.e., the extrapo¬ 
lated value of u t at which the saturated sand transport rate 
( Q s ) vanishes). This finding is consistent with the analytical 
model by Pahtz et al. [2013, 2014], supporting the hypothe¬ 
sis that the aforementioned four potentially most important 
relaxation mechanisms are all similarly relevant. 

This paper is organized as follows. First, we present 
the modeling technique we used to simulate aeolian particle 
transport in Section 2. Afterwards we show our numerical 
results in Section 3, which are then discussed and compared 
with the analytical model for the saturation length recently 
proposed by Pahtz et al. [2013, 2014] in Section 4. Finally, 
we draw conclusions in Section 5. 

2. Modeling technique 

In this section, we briefly describe the three-dimensional 
numerical model which we used to model sand transport. 
A more detailed description can be found in Carneiro et al. 
[2013], particularly its supplementary material. 

For the computation of the mean horizontal wind veloc¬ 
ity ( u ), the model uses the mixing-length approximation of 
the Reynolds-averaged Navier-Stokes equations (neglecting 
viscous wind shear), where the mixing length is given by 
k(z — h) with k = 0.4 being the von Karman constant and 
z — h characterizing the vertical distance from the top of the 
sand bed (z = h). This reads [Carneiro et al., 2013] 


OO 



The term containing /, the horizontal drag force per unit 
volume applied by the wind on the particles (the drag law 
by Cheng [1997] is used), takes into account that the wind 
speed is reduced due to continuous transfer of momentum 
from wind to particles. The integration of Eq. (1) starts 
at height z = h + z 0 , where z 0 = d/30 is the aerodynamic 
roughness of the sand bed as it would be in the absence 
of sand transport, corresponding to aerodynamically rough 
flow [Bagnold, 1941], However, during sand transport, the 
reduced wind speed at the top of the saltation layer corre¬ 
sponds to an increased roughness value (z*). 

It is important to note that Eq. (1) is applied to calcu¬ 
late the wind velocity profile at every time step because we 
assume that the flow adapts to local drag decelerations of 
the wind speed within the integration time (At = 0.005 s). 
This standard assumption led to several previous numeri¬ 
cal results in agreement with experiments [Carneiro et al, 
2011; Duran et al, 2012; Carneiro et al., 2013; Duran et al., 
2014b; Pahtz et al., 2015]. Nevertheless, we argue why it is 
reasonable to make such assumption in the following. 

There are actually two time scales involved. First, the 
time needed to transmit the drag force between fluid and 
particles. As the drag force is transmitted via collisions be¬ 
tween air molecules (which are extremely small) and the 


sand grains, this time scale is much smaller than the inte¬ 
gration time. The second time scale is related to the propa¬ 
gation of this perturbation to the entire system. As pertur¬ 
bations of a fluid typically travel at the speed of sound (in 
air, c « 340 m/s) and the maximum distance they need to 
travel to reach all relevant locations of the simulated salta¬ 
tion layer is of the order of lOOd = 0.02 m (the height of 
the saltation layer), the maximal time needed for this per¬ 
turbation to influence all relevant locations of the simulated 
saltation layer is around 0.0006 s, which is around a factor 
10 smaller than At and around a factor 10 3 smaller than the 
saturation time. Even if the necessary time for the flow to 
accommodate to a perturbation is larger than the necessary 
time a perturbation needs to travel to reach all relevant lo¬ 
cations, it is hard to imagine that this time comes any close 
to the saturation time. Moreover, since particle and flow 
velocity are of the same order of magnitude, also the length 
needed for the flow to adapt to the perturbation should be 
much smaller than the saturation length. 

Trajectories and velocities of particles are obtained from 
solving Newton’s equations of motion through the velocity- 
Stormer-Verlet scheme [Griebel et al., 2007], considering 
gravity and wind drag as the external forces acting on the 
particles. Interparticle contacts are modeled through a dis¬ 
sipative spring dashpot potential (coefficient of restitution, 
e = 0.65), while frictional contacts are neglected (no parti¬ 
cle rotation). The system dimensions are length x height x 
width = 50d x 400d x 7.5d, and 1410 particles with normally 
distributed diameters ( d p = (l±0.1)d) are simulated. Most 
of these particle constitute a bed of around twelve particle 
layers. This is sufficiently thick to suppress the reflection 
of shock waves from the dissipative (e = 0.5) bottom wall 
[Rioual et al., 2000, 2003]. The simulation top is open and 
the side boundaries periodic. In fact, particles never reach 
the top of the system. 

3. Results 

Using the model described in Section 2, we carried out 
simulations for typical sand transport conditions on Earth 



Figure 1. Time evolution of M/M s (red, dashed line), 
V/V s (green, dash-dotted line), and Q/Q s (blue, solid 
line) obtained from simulations of typical sand trans¬ 
port conditions on Earth (g = 9.81 m/s 2 , d = 200 g m, 
p p = 2650 kg/m 3 , pf = 1.174 kg/m 3 , v = 1.59 x 10 -5 
m 2 /s) with Ut = 1.2 m/s. The blue crosses correspond to 
the best fit of Eq. (5) to the data [ Q/Qs](t ) near satura¬ 
tion (|1 — Q/Qs\ < 0.2). The analogous plots for different 
values of u, look similar. 
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(g = 9.81 m/s 2 , d = 200 /im, p v = 2650 kg/m 3 , p / = 1.174 
kg/m 3 , v = 1.59 x 10~ J m 2 /s). For these conditions, we 
varied u* between two and nearly ten times the threshold 
shear velocity (ut = 0.195 m/s). For each u,, 15 runs were 
performed, starting from different initial condition. Indeed, 
while the sand bed in all samples was exactly the same at 
the start of each simulation (t = 0), ten particles with ve¬ 
locity v x = v z = 1 m/s were randomly placed sufficiently 
high above the surface (z > h + 20d). Each of these samples 
evolved in time towards the saturated state. The supple¬ 
mentary online material contains a movie showing the time 
evolution of a given sample (w» = 0.8 m/s) between 1 = 0 
and t = 0.35 s. 

From averaging the particle locations and velocities over 
the 15 samples corresponding to each u * and further over 
the horizontal ( x ) and lateral ( y ) direction, we obtain the 
vertical profiles of the local mass density (p) and the mass- 
weighted average particle velocity (v) the particle velocity at 
each time. The time evolution of M = pdz (red, dashed 
line), Q = p{v x )dz (blue, solid line), and V = Q/M 
(green, dash-dotted line) further obtained from these pro¬ 
files relative to their saturated values is plotted in Fig. 1 
for u t = 1.2 m/s (for different u», it looks similar). It can 
be seen that the time transient behavior of Q is similar to 
that determined in older studies [Spies et al. , 2000; Ma and 
Zheng, 2011]. Furthermore, one immediately recognizes that 
M relaxes significantly slower towards M s than V towards 
V s (this is also true for our simulations with other values 
of u t ). This means that our simulations do not confirm the 
hypothesis that drag is the dominant mechanism controlling 
sand transport saturation, which would instead require that 
M always relaxes much faster towards M s than V towards 
V s [Pahtz et al., 2014]. Moreover, this observation can be 
used to extract the saturation length (L s ) from the time 
evolution of Q, as we explain in the following. 

Mathematically, L s has its origin in the mass conserva¬ 
tion equation, which in its local form reads 



u„/u t 


Figure 2. The saturation length ( L s ) as a function 
of Ut /ut (symbols) obtained from simulations of typi¬ 
cal sand transport conditions on Earth (g = 9.81 m/s 2 , 
d = 200 pm, p„ = 2650 kg/m 3 , p/ = 1.174 kg/m 3 , 
v = 1.59 x 10 _ m 2 /s), whereby the dashed line cor¬ 
responds to L a = 2.30ut/g — 0.13 m and the dotted line 
to L a = 1.76 Ut/g. The inset shows L s as a function of 
V^/g for the same conditions, whereby the solid line cor¬ 
responds to L s = 0.481 ^s/g. The error bars correspond 
to the 95%-confidence intervals obtained from the best 
fits to Eq. (5). 


solution 



(5) 


dp dp(v x ) dp{v y ) dp(v z ) _ 
dt + dx + dy + dz 

Indeed, L s is typically defined through a first-order Taylor 
expansion of the rate of relaxation (T(Q)) of Q around Q 3 
[Andreotti et al, 2010; Pahtz et al, 2013, 2014], 


d Q 
dx 


T(Q) 


Qs-Q 

L s 


(3) 


which, using T = [p(v z )]{h), is the height integration 
(f//° -dci) of Eq. (2) for steady and laterally homoge¬ 
neous conditions ( d/dt = d/dy = 0). Since r(Q s ) = 0, 
L 3 corresponds to the negative inverse Taylor coefficient 
( — [T^Qs)] 3 ). By definition Eq. (3), describing the spatial 
relaxation of Q towards Q s , is only applicable near satura¬ 
tion (|1 — Q/Qs | 1). 

Since our simulations correspond to spatially and later¬ 
ally homogeneous conditions ( d/dx = d/dy = 0), height 
integration of Eq. (2) yields 


d Q ~ dM_ _ ^ Qs-Q 

Vdt ~ dt W> ~ L s 


(4) 


where the approximation on the left hand side uses that M 
relaxes significantly slower towards M s than V towards V s . 
From comparison between Eqs. (3) and (4), it is apparent 
that the saturation in t in our simulations is equivalent to 
a saturation in * if x = f 0 Vdt' (i.e., d/dx = Vd/dt) is 
used to relate them. In fact, fitting (nonlinear least squares 
method) Q s , x 0 , and L s to best agreement with the analytic 


of Eq. (4) (since Q s does not depend on t for constant wind 
shear) allows us to determine L s from our simulations (see 
blue crosses in Fig. 1). 

Fig. 2 shows L a as a function of Ut/ut obtained from our 
simulations, whereby the error bars correspond to the 95%- 
confidence intervals obtained from the best fits to Eq. (5). 
It can be seen that L a remains nearly constant between 2 ut 
and 4 ut, qualitatively consistent with measurements [An¬ 
dreotti et al, 2010]. However, L s increases with u* when 
Ut > dut- This increase within the error bars follows the 
scaling relation 

L 3 = aV s 2 /g, (6) 

as shown in the inset of Fig. 2, where a = 0.48. Interest¬ 
ingly, while the scaling L s oc V 2 /g describes the data with 
Ut > dut very well, the scaling L s oc ut/g does not (see 
the dotted line in Fig. 2, corresponding to the best fit of 
L s oc u 2 /g to data with Ut/ut > 4), indicating that V s and 
not Ut is the relevant parameter controlling L s . Only when 
one allows a small offset, a good fit can be obtained. This 
is shown by the dashed line in Fig. 2, which corresponds 
to L a + L so oc u 2 /g, where the offset ( L so = 0.13 m) is ex¬ 
pected to have a complex dependency on particle and wind 
parameters (except u»). 

4. Discussion 

In this section, we first briefly describe the analytical 
model for the saturation length proposed by Pahtz et al. 
[2013, 2014] in Section 4.1. We then compare the model 
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predictions with our numerical results shown in Fig. 2 and 
with the scaling in Eq. (6). 


4.1. Analytical model by Pahtz et al. [2013, 2014] 

The analytical model for the saturation length by Pahtz 
et al. [2013, 2014] takes into account that M and V relax 
towards M s and V 3 . respectively, due to different relaxation 
mechanisms. Changes in M are controlled by the ejection 
of bed particles in particle-bed collisions [Kok et al., 2012], 
while changes in V are driven by the acceleration of trans¬ 
ported particles due to fluid drag and their deceleration in 
particle-bed collisions. On the one hand, the changes of V 
were explicitly modeled within the momentum balance: the 
fluid drag acceleration through the fluid drag law by Julien 
[1995] for natural sand and the deceleration in particle-bed 
collisions by means of a Coulomb friction law, assuming 
a proportionality between average horizontal and vertical 
forces acting on transported particles. This assumption is 
well established in the literature as it leads to the experi¬ 
mentally [ Creyssels et al, 2009] and numerically [Kok and 
Renno, 2009; Pahtz et al., 2012] confirmed relation 

M s = ^(u 2 ,-u 2 ), (7) 


where p m 1 (from experimental data [Creyssels et al, 2009; 
Pahtz et al., 2012]) is the associated Coulomb friction coeffi¬ 
cient. On the other hand, changes in M were not explicitly 
modeled, but implicitly accounted for in the parameter 


V, d M 

CM = W s W {Va) ’ 


( 8 ) 


which describes the relative change of M with V near the 
saturated regime [Pahtz et al, 2014]. Moreover, taking into 
account that not only M and V, but also the average wind 
speed (U) relaxes towards its saturated value, led to the in¬ 
troduction of another parameter cxj, describing the relative 
change of U with the bed fluid shear velocity (i.e., the value 
of the fluid shear velocity at the bed, which is smaller than 
w» due to momentum transfer from fluid to particles). The 
parameters cm and cu were by far the most uncertain model 
parameters as they were the only ones not determined by 
measurements, but instead by theoretical arguments (which 
led to cm ~ cu ~ 1 for aeolian sand transport) [Pahtz et al, 
2013, 2014]. However, in the limit of large fluid shear ve¬ 
locities (u*/ut 1), the model becomes independent of cu ■ 
Indeed, in this limit, the final model equation for L s reads 
[Pahtz et al, 2014], 

Ls= (2±cm)c 2L VI 
cmr g 

where c v m 1.3 (from experimental data [Rasmussen 
and S0rensen, 2008; Creyssels et al., 2009; Pahtz et al., 
2013, 2014]) is the saturated value of p(v 2 )dz / (MV 2 ). 


4.2. Comparison between analytical and numerical 
model predictions 

It can be seen that Eq. (2) predicts L s oc V 2 /g, which re¬ 
sembles the scaling in Eq. (6) obtained from our simulations. 
Since both equations are only valid for sufficiently large val¬ 
ues u t /ut, this resemblance supports the analysis by Pahtz 
et al. [2013, 2014], In this analysis, the four potentially most 
important relaxation mechanisms are all accounted for (ejec¬ 
tion of bed particles and particle deceleration in particle-bed 
collisions, fluid drag acceleration of particles, relaxation of 


the fluid speed), as described in Section 4.1. Since neglecting 
any of them entirely changes the model predictions [Pahtz 
et al, 2014], the resemblance between the analytical and 
numerical model predictions suggests that these four relax¬ 
ation mechanisms are all similarly relevant. This is just 
another indication that Ld is not the dominant length scale 
controlling sand transport saturation. 

However, one must also note that there are differences be¬ 
tween these models. First, the qualitative model predictions 
for small it»/ut slightly differ from each other. While the 
numerical model predicts that L s remains nearly constant 
between 2 ut and 4itt, the analytical model predicts a slight 
increase with it* [Pahtz et al., 2013, 2014]. This might be 
a result of the aforementioned uncertainty of the parameter 
cu- Second, while the analytical model predictions are con¬ 
sistent with the measurements by Andreotti et al. [2010], the 
numerical model predictions shown in Fig. 2 are not. This 
can be entirely linked to differences in the model parameters 
p and cm, as we explain in the following. First, from fitting 
Eq. (7) to our numerical data (it fits very well, not shown), 
we obtain p m 2, in contrast to p m 1, which Pahtz et al. 
[2013, 2014] obtained from experimental data. Second, since 
M relaxes much slower towards M s than V towards V s in 
the simulations (see Fig. 1), cm becomes very large, while 
cm ~ 1 was estimated by Pahtz et al. [2013, 2014] from 
theoretical arguments. In fact, using pm 2 and cm -4 oo 
in Eq. (2) yields the prefactor (2 + cm)c v /(cmp) 0.65 
close to the prefactor a = 0.48 in Eq. (6). This means the 
numerical model and the analytical model seem quantita¬ 
tively consistent with each other since the differences in the 
parameters p and cm are likely the results of simplifications 
in the numerical model. For instance, the difference in the 
value of p between simulations and measurements can be 
linked to the interparticle contact model, which is known 
to have considerable influence on the frictional behavior of 
solids [Campbell, 2006]. Indeed, the model neglects particle 
rotation and uses rather soft particles (stiffness k = lkg/s 2 ), 
which allows particle overlaps of about 20%, while in reality 
the stiffness is several orders of magnitude larger. Also, the 
coefficient of restitution used in our simulations (e = 0.65) 
might have been too small. Comparable studies usually use 
larger values [e.g., Duran et al, 2012, e = 0.9] and obtain a 
Coulomb friction coefficients near unity. 

5. Conclusion 

We simulated aeolian sand transport using DEM simu¬ 
lations. From these simulations, we obtained the satura¬ 
tion curves in Fig. 1 of the total mass of particles trans¬ 
ported above the sand bed (M), their average velocity (V), 
and the associated sand flux (Q = MV). These numerical 
data indicate that M saturates much slower than V, chal¬ 
lenging the widely-accepted hypothesis that the drag length 
(Ld = ( Pp/pf)d ) is the dominant length scale controlling ae¬ 
olian sand transport saturation, which would require the op¬ 
posite, namely that M saturates much faster than V. Since 
Ld does not change with u*, the same hypothesis is fur¬ 
ther challenged by the numerical data in Fig. 2 showing 
that the saturation length ( L a ) significantly increases with 
the wind shear velocity (u*) for medium and strong winds 
(«* > 4 ut). Moreover, this increase follows the scaling rela¬ 
tion L s oc V 2 /g (see inset of Fig. 2), qualitatively consistent 
with the limit u*/ut 1 of the recently proposed analytical 
model by Pahtz et al. [2013, 2014]. In Section 4, we showed 
that this analytical model also predicts the proportionality 
factor in L s oc V 2 / g to be about 0.65, which is close to 
the numerically obtained value a = 0.48. This adds another 
piece of doubt on a dominating role of Ld because the model 
accounts for the four potentially most important relaxation 
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mechanisms (ejection of bed particles and particle decelera¬ 
tion in particle-bed collisions, fluid drag acceleration of par¬ 
ticles, relaxation of the fluid speed), and neglecting any of 
them entirely changes the model predictions [Pdhtz et al ., 
2014], 

The scaling relation L s oc V a /g, found for medium and 
strong winds, might itself become an important step to¬ 
wards modeling of dune and dune field evolutions in sand 
storms. For this purpose, one would need a model predict¬ 
ing V s . In fact, there are a considerable number of analytical 
models predicting Vs as function of u » [e.g., Bagnold, 1941; 
Kawarrmra, 1951; Owen, 1964; Bagnold, 1973; Kind, 1976; 
Lettau and Lettau, 1978; Ungar and Haff, 1987; Sgrensen, 
1991, 2003; Duran et al, 2011; Pdhtz et al, 2012; Lammel 
et al, 2012], some of them might be applicable to sand storm 
conditions. 

Finally, it is worth to note that aeolian dunes are often 
superimposed by ripples. Compared to the flat sand bed 
condition present in our simulation, the presence of such 
ripples leads to a strong increase of the aerodynamic rough¬ 
ness (z*), which corresponds to a smaller wind and thus 
saturated particle velocity (14) in the saltation layer. The 
relation L a oc V a /g, found for medium and strong winds, 
would thus imply that the formation of superimposed rip¬ 
ples on the surface of aeolian dunes is associated with a 
simultaneous decrease of L s . 
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